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ABSTRACT 

Over the last decade there has been immense progress in the follow-up of short and long gamma-ray bursts (GRBs), resulting in a 
significant rise in the detection rate of X-ray and optical afterglows, in the determination of GRB redshifts, and of the identification of 
the underlying host galaxies. Nevertheless, our theoretical understanding of the progenitors and central engines powering these vast 
explosions is lagging behind, and a newly identified class of ultra-long GRBs has fuelled speculation on the existence of a new channel 
of GRB formation. In this paper we present high signal-to-noise X-shooter observations of the host galaxy of GRB 130925A, which 
is the fourth unambiguously identified ultra-long GRB, with prompt 7 -ray emission detected for ~ 20 ks. The GRB line of sight 
was close to the host galaxy nucleus, and our spectroscopic observations cover this region along the bulge/disk of the galaxy, and a 
bright star-forming region within the outskirts of the galaxy. From our broad wavelength coverage, we obtain accurate metallicity and 
dust-extinction measurements at the galaxy nucleus and at an outer star-forming region, and measure a super-solar metallicity at both 
locations, placing this galaxy within the 10 — 20% most metal-rich GRB host galaxies. Such a high metal enrichment has significant 
implications on the progenitor models of both long and ultra-long GRBs, although the edge-on orientation of the host galaxy does not 
allow us to rule out a large metallicity variation along our line of sight. The spatially resolved spectroscopic observations presented 
in this paper offer important insight into variations in the metal and dust abundance within GRB host galaxies. However, they also 
illustrate the need for integral field unit observations on a larger sample of GRB host galaxies of a variety of metallicities to provide 
a more quantitative view on the relation between the GRB circumburst environment and the galaxy-whole properties. 

Key words, gamma-ray burst: individual: GRB 130925A - galaxies: individual: GRB 130925A - galaxies: abundances - dust, 
extinction 


1. Introduction 


The collapsar model (Woosley 1993[ Woosley & MacFadyen| 
1999[ [Hirschi et aTj 2005 [ ’^on & Langer| 2005 [ Woosley 


& Heger 


2006|l has been largely successful at describing the 


overall observed properties of long-duration gamma-ray bursts 
(GRBs), and has been validated by the detection of supernovae 
coincident with the GRB explosions* (GRB 980425/SN1998bw; 
Galama et al., 1998, GRB030329/SN2003dh; Hjorth et al., 2003, 
Stanek et al., 2003). Nevertheless, the model in its simplest 
form has faced challenges; notably the recent rise in solar or 
super-solar metallicities measured in the environments of GRBs, 
which is in contention with the metallicity upper lim it of Z ^ 
0.3 Zq imposed by the standard collapsar models (Woosley 


1993 1 MacFadyen & Woosley 1999[ Woosley & MacFadyen 


* Based on observations taken at ESOWLT, Programme ID 09 LA- 
0703 and 092.A-023I 


1999 


2006 


Hirschi et al.| 2005[ [Yoon & Langer[ 2005 [ Yoon et al. 


Woosley & Heger[ |2006|l, and more recently, the detec 


tion of a possible new class of ultra-long duration GRB, which 
requires a longer-lasting central engine than can be provided by 
the Wolf-Rayet (WR) progenitors believed to produce standard 
long GRBs. 

There have thus far been four clear examples of ultra-long 
duration GRBs, all with continuous prompt 7 -ray emission last¬ 
ing from thousands to tens of thousands of seconds rather 
than the more typical 10 — 100 s observed in standard long- 

[20TT]i, 


duration GRBs. These are GRB 101225A (Thone et al. 


GRB 111209A dStratta et al.l|20131|Gendre et al.l|20131|Greiner| 
[eTaTl [20T5] |, GRB 121027A 1Hou et al.| |2014) , and more re¬ 
cently, GRB 130925A, which had the longest duration prompt 
emission thus far detected in a GRB (~ 20 ks; |Bellm et al.| 


* For the remainder of the paper we will use the term GRB to refer to 
the long-duration class of events unless specified otherwise. 
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|MT4l|Evans et al.||20T4l[Piro et al.[|20T4l|Zhao & Sha^|20T4l l. 

The difficulty in accounting for the source of the long-lasting 7- 
ray emission, as well as notable differences in the X-ray light 
curve properties, have led to the suggestion that these ultra-long 
GRBs are a new class of transient, with progenitors that are dis¬ 
tinct from those that produce standard long GRBs (Gendre et aH 


predications of sub-solar progenitors. However, there is thus far 
limited evidence to suggest that GRBs reside within metal-poor 


regions of their overall metal-rich galaxies (although see Per- 


|2013[[Boer et al.||2015[ [Levan et al.l|2014t [Greiner et al.||2015| r. 

However, other GRBs with evidence of long-lasting central en¬ 
gine activity, such as short periods of emission spread over long 
timescales (~ 1000s), or long-lasting flaring at X-ray energies, 
may imply that there is a continuous distribution in the prompt 
emission duration of long GRBs, with the ultra-long class rep¬ 
resenting the extreme end of the distribution (e.g.|Virgili et al.| 
20T3|[Zhang~erai:i[20T4l l. 

When considering those GRBs with continuous prompt 
emission for > lO'^s, which we will refer to as the ultra-long 
class of events, the main progenitors explored have been tidal 
disruption events (TDE), magnetars, and low-metallicity blue 

[20T^ 


|ley & Perleyl |2013| l. Some spatially resolved analysis of GRB 
host galaxies has been done to map out the optical stellar light 
([Eruchter et al.| |2006[ [Hammer et al.| 2006[ Svensson et al.j 
|2010j l and dust emission (Michalowski et al. |2014| l, providing 
some measure of the variation in environmental properties at the 
GRB site relative to the rest of the host galaxy. In a few cases 
spatially resolved spectral analysis has been possible by either 
using well-positioned single slit observations ([Hammer et al. 


2006| TTione et al.j 2008| [Levesque et al. 2011 , or in rarer cases. 


tntegral field unit (lEU) data ( [Christensen et al.[ [2008 [ [Thone et] 

[nij [2014| l. However, in most cases these studies have involved 
relatively metal-poor host galaxies. The one exceptio n was the 
super-solar host galaxy of GRB 020819B^, for which Levesque 


et al. (2010b[l attained spectra at the GRB explosion site and the 


supergiants (BSG) (Gendre et al. 


However, thus far the sample size of these events is too small 
to be able to strongly favour a single channel of formation. 
The unusually high host galaxy dust extinction measured along 
the line of sight to the most recently detected ultra-long GRB, 
GRB 130925A, introduces new diversity in the environmental 
properties associated with long, and especially ultra-long GRBs. 
With the typically rapid dissemination of sub-arcsecond GRB 
positions provided by the GRB-dedicated Swift mission ([Gehrels[ 


host galaxy nucleus, but found the metallicity to be the same in 
Levan et ah] [2014[ l. both regions. 

In this paper, we report on our UV through to NIR X- 
shooter spectroscopic observations ( [D’Odorico et al.[ 2006[[\fer^ 
net et al.| [201 l| l of the host galaxy of GRB 130925A, target- 


et al. 20041, launched in 2004, and the commissioning of (semi- 
) robotic telescopes equipped with near-infrared (NIR) facilities, 
such as the GRB optical and near-infrared detector (GROND; 
Greiner et al.[ 2008[l, the P eters automated IR imaging telescope 
(PAIRITEL; [Bloorn et al.[ [2006[) and the reionization and tran¬ 


sients IR camera (RATIR; Butler et al.[ [2012[ l, the detection 
of heavily dust-extinguished GRBs and dust-rich host galaxies 
has become more common. GRB 130925A was at a redshift of 
z = 0.347 ( [Sudilovsky et al.[|2013| l and had a host galaxy visual 
extinction of Ay= 5.0 ± 0.7 mag as measured from the after¬ 
glow spectral energy distribution (SED) (Greiner et al. 2014[l, 


which is one of the largest host galaxy visual extinctions mea¬ 
sured along the line of sight to any GRB, and certainly the largest 
along the line of sight to an ultra-long duration GRB. 

The host galaxies of heavily dust extinguished GRBs (i.e. 
Ay> 1 mag) typically have larger stellar masses, luminosities 
and dust masses than the hosts of optically bright GRBs ([Perley 


ing precisely to explore further the local environments of heav¬ 
ily dust-extinguished GRBs. Observations began 3.5 h after the 
Swift GRB trigger, and although by this time the heavily dust- 
extinguished afterglow was barely detected, our X-shooter slit 
was positioned such that it covered two spatially distinct star¬ 
forming regions within the host, one of which was coincident 
with the GRB location. 

The structure of the paper is as follows. In sectionj^we de¬ 
scribe briefly our GROND and X-shooter observations, and the 
X-shooter data reduction is described in section We present 
our X-shooter spectral data analysis and results in section]^ and 
in section]^ we discuss the implications that these results have 
on our understanding of the environments of ultra-long duration 
GRBs, and their progenitors. Throughout the paper, we assume 
a Om = 0.27, Ha = 0.73 cosmology, with Hubble constant 
Hq = 70 km s“^ Mpc“^, and all reported errors are Icr. 


2. GRB and host galaxy observations 

GRB 130925A was detected by numerous gamma-ray detectors 


et al.[ 2009 1 Kriihler et al.[ 201 1[ [Rossi et al.[ 2012[ Perley et 


al. 20131 Hunt et aL[ 2014[ Schady et al.| 2014|i, and there are 

now an appreciable number of GRB host galaxies with near-solar 

or super-solar metallicities (Prochaska et al. 

2007[ Levesque et 

al. 2010b[ Savaglio et al. 

2012| Elliott et a 

.[ 20131 Graham & 

Eruchter[ 2013[l. Althoug 

1 these more massive and metal-rich 


(|Fitzpatrick| 

2013 

Hurley et al. 

2013 Jenke| 2013 Savchenko 

[et al.[ [20131 

Suzu 

ki et al.[[2013 

1 , and the Swift burst alert tele- 

scope trigger (BA 

Barthelmy et al.[ 2005[l led to an immediate 


slew and subsequent X-ray afterglow detection with the X-ray 
telescope (XRT; [Burrows et al.| [2005 [) at the refined position 
RA(J2000.0)=02:44:42.91, Dec (J2000.0) = -26:09:10.8 with 


an error radius of 1.6 (Evans et al. 2013 1 . No optical afterglow 


GRB host galaxies do not represent the majority of the GRB 
host galaxy sample, they do call for certain details of the classi¬ 
cally accepted collapsar model to be reviewed, such as the effect 


of differential rotation (Georgy et al.[[2012|l, and binary system 

models ([Podsiadlowski et al.[ 2004[ Izzard et al. 

|2004[|Petrovic| 

[et al.[[2005 

van den Heuvel & Portegies Zwart 

20131, or alter- 


redward of the i?-band was detected, with the ultraviolet and op¬ 
tical telescope observations (UVOT; Roming et al. 2005| l taken 
just 157s after the GRB trigger, resulting in a u-band 3cr upper 
limit of u > 20.0 ( [Holland &Lien[[201 3[ l. Rapid-response obser¬ 
vations with the simultaneous 7 channel imager, GROND, began 
just 7 minutes after the GRB trigger (Sudilovsky et al. 2013[l, 


natively, for a re-assessment of how representative the galaxy 
whole, or even the GRB line of sight properties are of the condi 
tions within the GRB local environment. 

Erom an observational perspective, sharp galaxy metallicity 
gradients, or low-metallicity ‘bubbles’ within the GRB vicinity 
have been proposed (e.g. Campisi et al. 2009 [ 


and a decaying optical source was detected within t he X-ray af- 
terglow error circle in the i',z',J,H and K bands (Sudilovsky 
et al. 2013[ Greiner et al. 2()14| l. A detailed analysis of the XRT 


20111 Kocevski et 


al. 20091, in order to consolidate the more chemically evolved 


and dust-rich GRB environments measured in some cases with 


^ Although this GRB was originally designated GRB 020819A, it is 
in fact the second GRB to have been detected that day, and is thus 
020819B. GRB 020819A was detected by IPN satellites ~ 7 h prior 
to GRB 020819B. See http://www.ssl.berkeley.edu/ipn3/masterli.txt for 
details. 
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and GROND afterglow observations are given in [Evans et al.| 
( 2014| l and Greiner et al. ( 2014| l, respectively. 

Rapid response mode spectroscopic observations of 
GRB 130925A with the VLT/UVES spectrograph began 
50 minutes after the Swift trigger, and revealed faint continuum 
emission and two emission lines identified as Ha and [N ll] at a 
redshift of z = 0.347 ( [Vreeswijk et al.||2013| l. Owing to the lack 
of a bright afterglow, and thus of absorption lines, the UVES 
observations were aborted, enabling us to trigger our X-shooter 
ToO programme (ID 091.A-0703; PI: P. Schady). 

Our X-shooter spectroscopic observations began on 25-09- 
2013 at 07:43 UT, 3.5 h after the BAT trigger. Because of the 
disabling of the X-shooter atmospheric dispersion correctors in 
August 2012, the slit position was orientated to the parallactic 
angle in order to minimise atmospheric dispersion losses. The 
observations consisted of four nodded exposures in the sequence 
ABBA with exposure times of ~ 1500s taken simultaneously 
in X-shooter’s untraviolet/blue (UVB), visible (VIS) and near- 
infrared (NIR) arms. The average airmass during our observa¬ 
tions was 1.0 and the seeing was a median of 0.6". X-shooter 
spectroscopy was performed with slit-widths of l."0, 0."9 and 
0."9 in the UVB, VIS and NIR arm, respectively. The resolving 
power R = A/A A was determined from arc lamp calibration 
data taken with the same instrumental setup as our science ob¬ 
servations^. Prom this we get R ~ 5300 for the VIS and NIR 
arms, and R ~ 8800 for the VIS arm. These correspond to an 
instrumental velocity FWHM of ^ 55 km s“^ and ^ 35 km s“^, 
respectively. 

HST observations taken 20 days after the GRB indicate that 
the GRB lies close to the plane of the disk of what appears to 
be an edge-on spiral, just ~ 0."12 from the nucleus (~ 600 pc 
in projection) ( [Tanvir et al.||2013] l. A detailed discussion on the 
HST observations of the host galaxy will be reported in Tanvir et 
al. (in prep.). We also observed the field of GRB 130925A with 
the VLT/HAWK-I instrument in the Y, J, H and Ks bands 18 
days after the GRB trigger (ID 092.A-0231; PI: T. Kriihler), and 
the measured photometry was consistent with the source mag¬ 
nitudes measured seven days earlier with GROND, suggesting 
that these late-time GROND and HAWK-I observations were 
dominated by emission from the host galaxy. From our more 
sensitive HAWK-I data we measure host galaxy AB smagni- 
tudes of Jab = 20.70 ± 0.07 mag, Hab = 20.47 ± 0.07 mag 
and Kab — 20.20 ± 0.08 mag (Greiner et al. |2014| l, cor¬ 
responding to an absolute (K-corrected) iT-band magnitude of 
niHiAB) = -20.6 ± 0.08 mag. The HST F814W image with 
the position of the X-shooter slit overlaid is shown in Fig. [T] 


3. X-Shooter data reduction 

X-shooter data were reduced with the ESO/X-shooter pipeline 
v2.2.0 ( [Goldoni et al.j [2006| l, including the flat-fielding, order 
tracing, rectification and initial wavelength calibration from arc 
lamp data. During rectification, a dispersion of 0.2 A/pix and 
0.6 A/pix was used in the UVB, VIS and NIR arms, respectively, 
which minimises correlated noise while still maintaining suffi¬ 
cient spectral resolution to resolve lines down to ^ 25 km s“^ in 
the UVB and VIS arms (i.e. a velocity dispersion of 10 km s“^), 
and down to ^ 50 km s“^ in the NIR arm (i.e. a velocity disper¬ 
sion of 20 km s“^). Our own software was used for bad pixel and 
cosmic-ray rejection, for sky subtraction, and for frame shifting 
and adding (see [Kriihler et al.|2015[ for details). The X-shooter 

^ archive.eso.org/bin/qcl_cgi 



Fig. 1. HST/F814W image of the host galaxy of GRB 130925A. 
The image is 4"x4"centred on the GRB host galaxy. It is displayed 
with a linear greyscale ranging from 0 count s“^ pix“^ (white) to 
0.1 count s“^ pix“^ (black), and logarithmically spaced contours are 
shown in white. The position of the X-shooter’s VIS slit with width 
0."9 is indicated by the black dashed lines. The UVB and NIR slits 
have the same orientation. The position of the afterglow was determined 
from image subtraction between two HST/F814W observations taken 
20 and 48 days after the GRB, and is indicated by the red cross. It falls 
within the emission region SFl detected in our spectral data, which in 
the above image, corresponds roughly to the second dimmest contour 
plotted. The emission region, SF2, detected in our spectra and located 
l."2 from SFl is indicated at the SW corner of the outer contour plot 
by a blue ellipse. 


spectra were then flux calibrated with the nightly spectrophoto- 
metric standard star and optimally extracted within IRAF using 
a large aperture in order to include all the emission from the 
galaxy. The Icr error bars on each data point in our extracted 
ID spectra were derived via the error spectrum and propagated 
accordingly. 


4. Data analysis and results 

The X-shooter spectrum of GRB 130925A is dominated by 
emission lines from the host galaxy, with only a very weak con¬ 
tinuum detected, corresponding mostly to emission from the host 
galaxy, with a small contribution in the NIR arm from the af¬ 
terglow. This is verified by the continual GROND observations 
taken between 350s and 10®s after the GRB trigger, which show 
that at the epoch of our X-shooter observations, 3.5 h after the 
GRB trigger, the emission predominately originates from the 
host galaxy light, with the afterglow contributing only a small 
amount to the flux in the GROND J, H and K bands ( [Greiner ^ 
[ar||2014| . In both the 2D and the ID spectra, we clearly detect 
all the usual strong emission lines, and H/3, H 7 , [N It] and the 
[S ll](A6717 & A6731) doublet are all detected at better than 3 (t 
( see Fig. [^. In our 2D spectra, emission from at least two star¬ 
forming (SF) regions, SFl and SF2, is detected, corresponding 
to emission from the central disk/bulge of the galaxy, and from 
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Fig. 2. Balmer and forbidden emission lines from the host galaxy of GRB 130925A. In each figure, the top and middle panels show the X-shooter 
2D and ID spectra zoomed into the wavelength range of the line labelled in the top right comer of the middle panel. The bottom panel shows the 
residuals from our two-component fits to the ID spectrum. In the middle panel, we also show the best-fit Gaussians for each individual component 
(dashed grey lines) and composite model (red solid line). 


an HII region located ^ l."2 south-west of the galaxy nucleus 
(i.e. 6 kpc in projection; see Fig.[2l. The GRB afterglow position 
is coincident with the location of the weaker and extended SFl 
region. 

We extracted a single ID spectrum for each of the X-shooter 
arms using a large aperture of width 4."2, which contains the 
emission from both SFl and SF2. The emission profiles for most 
of the detected lines clearly deviate from a Gaussian, and we 
therefore fitted each emission line in our ID spectra with two 
Gaussians simultaneously, one per detected emission region. For 
line doublets we tied the Gaussian width of the lines correspond¬ 
ing to the same emission component, and in the case of the weak 
H 7 line, we fixed the Gaussian width and central line peak posi¬ 
tions for the two emission components to correspond to the best 
fit redshift and velocity line width determined from our fits to 
the Ha line. In our fits to the [O ll] and [S ll] doublets, and to 
the [O lll]A4959 line, we fixed the central line peak position of 
the weaker Gaussian component from SFl to correspond to the 
best fit redshift measured at SFl in our Ha line fits. 

Slit losses on our line emission flux measurements due to 
the finite width of the slit were estimated independently for each 
arm against the multi-band GROND host galaxy photometry, and 
were found to be consistent within errors. The corresponding 
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average slit loss correction was a multiplicative factor of 2.7 ± 
0.1. The spatial offset between the two detected emission regions 
introduces further uncertainty in our slit-loss correction, which is 
likely to not be the same for the two star-forming regions. Given 
the proximity of SF2 to the slit edge (see Fig. [^1, the respective 
slit losses might be notably larger than that at SFl. This will not 
affect greatly the dust extinction or metallicity that we measure 
at SFl and SF2, which is the focus of this paper. However it may 
impact our measured line fluxes (see Table [TJ and derived star- 
formation rate (SFR), such that they are over-estimated at SFl 
and under-estimated at SF2. 


We are unable to resolve and thus correct our measured 
Balmer line fluxes for Balmer absorption (e.g. |Savaglio et al.| 
|2005| ICovino et al.| |2006| ). However, for such a weak stellar 
continuum, we expect the effect of Balmer absorption on our 
line fluxes to be small, with other factors such as slit loss and 
dust extinction corrections being a more dominant source of un¬ 
certainty. Our best-fit slit loss corrected line fluxes are listed in 
Table [T] where the given errors include the uncertainties on the 
location of the Gaussian peak flux and the FWHM. Systematic 
photometric calibration errors and slit-loss correction errors are 
not included, which correspond to a further 5-8%. 
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Table 1. GRB 130925A host galaxy line fluxes at SFl and SF2 cor¬ 
rected for a Galactic reddening of E(B-V)=0.02 



Wavelength® 

Fa (SEl) 

Fa(SE2) 


A 

(10 erg cm ^ s ^) 

H 7 

4340.5 

0 4rj+u.:^fc) 

— 1 77 + 0.19 — 

' ' - 0.19 

H/3 

4861.3 

A 551 + 0-50 
4 - 0 - 1 -_o .44 

6.30+:“ 

Ha 

6562.8 

27.8077^); 

27.671+° 

[Oil] 

3726.0 

10.83ti-?l 

' -^^-0.10 

[On] 

3728.8 

12.55t7« 

11.84+:“ 

[O III] 

4958.9 


4 C7+O.12 
^•^^-0.11 

[O III] 

5006.8 

V 521 +0.11 

iz.yu_Q 

[Nil] 

6583.4 

c; Al + 0.30 
^•'J-‘-- 0.47 

'^•^'^-0.11 

[S II] 

6716.5 

O.UO_o,69 

0 40 + 0.05 

[S II] 

6730.9 

O.00 o ,45 

^•^^- 0.05 


Rest-frame vacuum wavelength 


4.1. Velocity dispersion and dynamical mass 


Our two-component spectral fits give best-fit Gaussian widths 
corresponding to velocity dispersions of 80 — 90 km s“^ at SFl, 
and 25 — 30 km s“^ at SF2, which is close to the velocity res¬ 
olution of our data. Such velocities are at the lower and higher 
^ 20 % of the distribution of velocity dispersions observed in 
GRB host galaxies at redshifts z < 1 (Kriihler et al.| 2015|l. 


From our best-fit emission line central wavelengths we get he¬ 
liocentric'^ redshifts of zsfi = 0.3479 and zsf 2 = 0.3483 
at SFl and SF2, respectively, corresponding to a separation in 
velocity space of 90 km s“^. Within the projected distance of 
^ 6 kpc of SF2 from the galaxy nucleus, a velocity dispersion 
of 90 km s“^ corresponds to a logarithmic dynamical mass of 
\og[Mdyn/MQ] ~ 10.1. 


4.2. Balmer decrement and SFR 


Assuming a temperature T=10'*^ K and a density of 10^ — 
10^ cm“^, we compare the observed Ha/H/3 Balmer decrement 
to the intrinsic, unextinguished ratio of 2.86 in order to measure 
the host galaxy visual extinction. Ay, at SFl and SF2. Assuming 
a Calzetti dust attenuation law ( |Calzetti et al.||2O0O| l, we measure 
Ay=2.4 ± 0.9 mag and Ay=1.5 ± 0.1 mag at SFl and SF2, re¬ 
spectively. Comparing the measured HalW'y ratio to the intrinsic 
ratio of 6.11 results in visual extinctions of Ay=0.7 ± 0.2 mag 
and 2.3 ± 0.6 mag at SFl and SF2, respectively, which are con¬ 
sistent at 2tT with the Ay values determined from the measured 
Ha/H/? ratio. Given the relative weakness of the H 7 emission 
line, and thus the larger uncertainty in the corresponding flux 
density measurement, we choose to use the Ay values derived 
from the Ha/H/3 ratio for SFl and SF2. This gives us a dust- 
corrected Ha luminosity of (7.1 ± 5.2) x lO'^^ erg s“^ and 
(3.6 ± 0.3) X lO'^^ erg s“^ at SFl and SF2, respectively, coiTe- 
sponding to SFRs of 3.4 ± 2.5 Mq yr“^ and 1.7 ± 0.1 Mq yr“^ 
based on the formalism described in |Kennicutt| d 199811, but con¬ 


verted to a Chabrier IMF ( |Chabrier7 2003|l, which re duces the 
SFR by a factor of ~ 1.7 (e.g. Forster Schreiber et al.||200^ . 


The heliocentric correction in the direction of GRB 130925A is 
10 km for our observations. 


4.3. Host galaxy metallicity 


There are a number of available empirical and theoretical metal¬ 
licity diagnostics, each with their own strengths and drawbacks. 
Owing to the range in emission line ratios and physical assump¬ 
tions (e.g. ionisation parameter) that go into the different gas- 
phase metallicity measurements, certain diagnostics are more 
appropriate for sources within a given metallicity or redshift in¬ 
terval. Nevertheless, which of the metallicity diagnostics is used 
is typically limited by the emission lines detected for the source 
in question, and thus it is not always possible to use the most ac¬ 
curate metallicity diagnostic. As a result of the very large wave¬ 
length range available with X-shooter, we measure all the rele¬ 
vant emission lines (Table required from the two star forming 
regions to apply several of the more accurate metallicity diag¬ 
nostics. 


In a detailed investigation on the accuracy and cross¬ 
calibration between the more popular theoretical and empirical 
metallicity diagnostics, |Kewley & Ellison] ( |2008[ ) derive metal¬ 
licity conversions between eight of the strong-line metallicity 
calibrations, allowing the metallicities to be converted to the 
same base calibration. Nevertheless, such conversions are based 
on the average difference between metallicity diagnostics, and 
are therefore designed to be used on galaxy samples, rather than 
on individual sources. Below we therefore derive the host galaxy 
metallicity at SFl and SF2 for a number of diagnostics, and give 
the metallicities using the native metallicity diagnostic calibra¬ 
tion, as well as the metallicities converted to the |Kobulnicky & 


Kewley p0()4|l scale, applying the conversions from Kewley & 


Ellison| (|2008|l. The latter conversion to a base scale facilitates 


a comparison between the derived metallicities and solar m etal- 
licity, which we take to be log(0/H) -f 12 = 8.69 ± 0.05 (Al- 
lende Prieto et al. 2001[ |Asplund et ah |2009 l. The R 23 diag- 
nostic from |Kobulnicky & Kewley| ( |2004[ ) iscalibrated against 
this same solar metallicity (see their footnote 2 ), which is why 
we choose this metallicity diagnostic as our base scale. 


[Kewley & Ellis ori] ( |2008|) found that the N 202 metallicity 
diagnostic as calibrated by [Kewley & Dopitaj ( [2002[ l was the 
most reliable for systems with log[N Il/O llj > —1.2 (corre¬ 
sponding to 12 + log[0/H] > 8.4) because of its weak depen¬ 
dence on ionisation parameter within this metallicity regime. Eor 
12 + log[0/H] < 8.4 the dependence on ionisation parameter in¬ 
creases, and thus in such a case [Kewley & Ellison| ( [200^ recom¬ 
mend using the R 23 method, which estimates the ionisation pa¬ 
rameter from the observed [O ni]/[0 ll] line ratio. The drawback 
of the i ?23 method is that it has a lower and an upper branch solu¬ 
tion. However, the degeneracy between the two solutions can be 
broken using other metallicity diagnostics, such as in this case by 
the N202 diagnostic, whereby the lower branch solution is ap¬ 
propriate for sources with a flux line ratio log[N Il/O ll] < —1.2. 


In the case of the host galaxy of GRB 130925A, we mea¬ 
sure log[N Il/O ll] = —0.7 at SEl and SE2, after correcting 
for host galaxy dust reddening. The [Kewley & Dopita] ([2002|l 


N202 diagnostic is thus appropriate, and gives us a metal¬ 
licity of 12 -F log[0/H]N202 = 8.77ti(;;;2 at SEl, and 12 -F 



reported errors include the uncertainties in the line flux mea¬ 
surements and host galaxy reddening correction derived from 
the Balmer decrement, but not the systematic calibration error 
of ~ 0.1 dex. 
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|Kewley & Ellison (|2008|l also found the 03N2 and N2 diag¬ 
nostics from Pettini & Pagel| ( [2004| l to be fairly robust, with the 
added advantage that these methods are fairly insensitive to the 
dust reddening correction applied as a result of the proximity in 
wavelength space of the lines used. This feature makes these two 
diagnostics favourable for high redshift (z > 1) sources, where 
there may be a large uncertainty in the measured Balmer decre¬ 
ment. Although this is not a significant issue in our data, we 
nevertheless calculate the metallicity measured with these two 
diagnostics, and find that 12 + log[0/H]o3N2 = 8.50 ± 0.03 and 
12-flog[0/H]N2 = 8.60±0.04atSFl, and 12-f log[0/H]o3N2 = 
8.3 5 ± 0.01 and 12 -f log [ 0/H]N 2 = 8.33 ± 0.01 at SF2. In 
the |Kobulnicky & Kewley| ( |2004[ ) scale, the metallicity at SFl 
is 12 -I- log[O^Jo 3 N 2 = 8.85 ± 0.03 and 12 -f log[0/H]N2 = 
8.90 ± 0.04, and it is 12 + log[0/H]o3N2 = 8.73 ± 0.01 and 
12-|-log[0/H]N2 = 8.72 ±0.01 atSF2. These values correspond 
to respective 03N2 and N2 metallicities of 1.4^0 and 1.6Zq 
at SFl, and to I.IZq at SF2 in both the 03N2 and N2 metal¬ 
licity diagnostics. When considering the respective systematic 
calibration errors of ^ 0.2 dex and ~ 0.4 dex for the 03N2 and 
N2 methods, the derived metallicities are consistent at Icr with 
those derived with the N202 calibrator. 

Given that the R 23 method is frequently used, for complete¬ 
ness we also use this diagnostic to measure the metallicity in 
SFl and SF2, applying the calibration from Kobulnicky & Kew- 


ley 


( 2004) 1, which simultaneously fits the ionisation parameter 
and gas-phase metallicity. This gives us an upper-branch solu¬ 
tion of 12 ± log[ 0 /H]i ^,3 = 8.70lJ^ ie ^ 0 ) at SFl, and 
12 ± log[ 0 /H]i ^33 = 8.90 ± 0.01 (1.6Zq) at SF2, which is fully 
consistent with the N202, 03N2 and N2 derived metallicities 
when systematic errors are taken into account. 
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Fig. 3. SED of the host galaxy of GRB 130925A and best-fit galaxy 
model (solid line) in the observer frame. Filled black circles represent 
photometric measurements taken with GROND in g'r'i'z' and HAWK- 
I in JHK. 


4.5. Combined versus individual spectral analysis 

We re-did the above analysis on the individual SFl and SF2 re¬ 
gions by using smaller apertures to extract the ID spectra, cen¬ 
tred on the peak flux of each emission region. The apertures we 
used were l."4 and 2."1 wide for the SFl and SF2 regions, re¬ 
spectively. From the line fluxes measured from the separate ID 
spectra, we derived Balmer decrement corrections and metallic¬ 
ities at SFl and SF2 fully consistent with the values measured 
from our combined spectral analysis, although the error bars 
corresponding to the metallicity at SFl increased by a factor of 
around three. 


A solar or above metallicity is also supported by the lack of 
emission detected at the location of the [O lll]A4363 emission 
line down to a 3 ct upper limit of < 1.4e“^^ erg cm“^ s, whereas 
we would expect to detect emission from this line within a low- 
metallicity environment. 


4.4. Galaxy spectral energy distribution 

Using late time g'r'i'z'JHK GROND and JHK VLT/HAWK- 

I data taken 11 days and 18 days after the GRB trigger, respec¬ 
tively ( |Greiner et al.| |2014| l, we created the host galaxy SED. 
There was no evidence of further fading in the JHK filters be¬ 
tween these two epochs, indicating that we were detecting the 
host galaxy at these times with a negligible contribution from 
the afterglow. The flux in the optical bands plateaued at 2 days 
after the GRB trigger, with GROND observations taken 2, 4 and 

II days after the trigger yielding consistent magnitudes ( |Greiner| 
et al. |2014|l. We modelled the host galaxy SED using the spec 


tral fitting package LePHARE ( |Arnouts et al. 1999) Ilbert et 
[ah] |2006| l, which is a population-synthesis-based fitting proce¬ 
dure. We used the |Bruzual & Chariot) ( 2003] l galaxy templates, 
which include emission lines and prescribed reddening and pa¬ 
rameters therein, and we assumed a )Chabrier) P003)l IME (see 
Fig. 0. We measure a host galaxy logarithmic stellar mass of 
log[M*/Mo] = 9.5 ± 0.2 and SER of 2.4^^;® Mq yr"! for 
a host galaxy reddening of E{B — V) = 0.2, assuming again 
a Calzetti dust attenuation law ( jCalzetti et al.) )2000) l. This is 
consistent within the errors with the Ha luminosity-derived SER 
(section [4.2) ). 


5. Discussion 

5.1. Dust and metal abundances in GRB host galaxies 

The host galaxy of GRB 130925A is remarkable in terms of the 
amount of dust attenuation detected along the disk of the galaxy 
(Ay= 2.4 ± 0.9 mag), as well as in the outer regions of the 
galaxy (Ay= 1.5 ± 0.1 mag). The host galaxy dust-extinction 
measured along the GRB line of sight from the afterglow SED 
(Ay=5.0 ± 0.7 mag) is the second largest thus far measured 
in a GRB afterglow®, not considering Ay lower limits. The 
relatively low redshift of this GRB makes the high afterglow 
extinction all the more unusual, given that more heavily dust- 
extinguished GRB afterglows, and correspondingly dust-rich 


GRB host galaxies, are more numerous at high redshift (Perley 
)et al.))20T3))Kruhler et al.))2015) ). Whereas at z < 1 around 20% 
of GRB host galaxies have a visual extinction Ay > 1.2 mag, 
and less than 10% have Ay > 1.5 mag, at redshifts z = 1 — 2, 
^ 60% of GRB host galaxies have Ay > 1.2 mag, and over 40% 
have Ay >1.5 mag ( )Kruhler et al. 2015) 1. Prior to the commis¬ 
sioning of NIR, rapid-response facilities, such as GROND and 
PAIRITEL, the afterglow of this GRB at wavelengths longward 
of X-rays would likely not have been detected. 

The proximity of the GRB to the nucleus of the galaxy, and 
more importantly, the seemingly edge-on orientation of the host 
galaxy implies that the GRB line-of-sight crossed a particularly 
dusty region of its host, which may, at least in part, account 
for the unusually high GRB afterglow visual extinction. Nev- 


® The afterglow of GRB 070306 currently has the largest me asured 
host galaxy, afterglow visual extinction, with Ay= mag (Jaun- 

sen et al.)|20081|Kruhler et al.)|201 lb. 
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ertheless, it is noteworthy that GRB 130925A is only the fourth 
GRB host galaxy at z < 0.5 to have solar or above metallic- 
ity, and that a high metallicity is measured within the disk, and 
in a star-forming region in the outskirts of the galaxy. These 
observations make the host galaxy of GRB 130925A only the 
second super-solar host galaxy where spatially resolved spectro¬ 
scopic observations have been possible; the first being the host 
galaxy of GRB 020819B at z = 0.411 ( [Levesque et al.] 2010b| l. 
Furthermore, unlike the spectroscopic observations of the host 
galaxy of GRB 020819B for example, the unique broad wave¬ 
length range of X-shooter covers all emission lines necessary to 
be able to measure the metallicity relatively accurately. Acquir¬ 
ing such measurements for metal- and dust-rich GRB host galax¬ 
ies is challenging as a result of the initial difficulty in detecting 
the afterglow along such highly dust extinguished lines of sight, 
and thus acquiring an arcsecond GRB position, and secondly, in 
then obtaining sufficiently high signal-to-noise host galaxy spec¬ 
tra to measure the Balmer decrement and galaxy metallicity with 
any degree of accuracy. 

There are several similarities in the properties of the host 
galaxies of GRB 020819B and GRB 130925A, including a com¬ 
parably large visual dust extinction at the galaxy nucleus of 
Av^= 1.5 — 2.0 mag as derived from the measured Balmer decre¬ 
ment, and a super-solar metallicity at the nucleus and outskirts 
of the galaxy. The stellar mass, SFR and metallicity of the host 
galaxy of GRB 020819B were, however, somewhat larger than 
that of host galaxy of GRB 130925A (i.e. log M*/M^ ~ 10.5, 


SFR^ 15Mq yr-\ log(0/H) + 12 ~ 9.0 in the Kewley & 


Dopita||2002 scale; [Savaglio et al. 2009| Kiipcii Yolda§ et al. 


flares that were initially believed to originate from a TDE ( Bur- 
|rows et |2013|). The X-ray spectrum was also unusually soft 


(/3x ~ 3) and an emission component other than the standard 
external forward shock emission is required to account for the 
X-ray emission (Bellm et al. 2014| Evans et al. 2014 Piro et al. 


|2014| l. [Evans et al. 1 2014| l and Zhao & Shao ( |2014| l both argue 
that the late-time X-ray emission is the result of dust-scattering 
of the prompt emission down to X-rays, whereas Bellm et al.j 
( 2014]l and|Piro et al.| ( |2014] ) find that, similar to GRB 101225A 


(Thone et al. 2011 although see Levan et al., 2013), a black- 


body component in addition to the typical power-law external 
shock emission component, provides a good fit to the data (see 
also Basak & Rao 2015 for a two blackbody plus power-law 
model). Despite differences in the interpretation of the X-ray 
afterglow emission mechanism, a general consensus is that the 
slight offset between the GRB position and the host galaxy nu¬ 
cleus, and more significantly, the duration of the X-ray emission, 
favour a collapsar rather than TDE. Eurthermore, both [Evans etj 


al. (2014 1 and Piro et al. ( 2014| l find that the low luminosity of 
the external shock emission implies a low-density circumburst 
environment, which would be consistent with low mass-loss ex¬ 
pected from a BSG. However, such a progenitor model would 
also predict a low-metallicity environment (Gendre et al. [2013 


Nakauchi et al.[[2OT3]l, and in the case of GRB 130925A, 


Piro et[ 


al. (2014 1 proposed a Pop III progenitor star. This is in contrast 


2010[[Levesque et al.[[2010b[[Svensson et al. 2010| l. An impor- 
tant observational distinction between the two galaxies is that 
the host galaxy of GRB 020819B is face-on, and the GRB oc- 
cutTed within the outskirts of the galaxy. Spectroscopic observa¬ 
tions of this latter host galaxy thus probe smaller radial length 
scales than is the case in our observations of the host galaxy of 
GRB 130925A. The uncertainty on how representative the spec¬ 
troscopic measurements are of the conditions nearby to the GRB 
are thus larger in the case of GRB 130925A ( jNiino et aLl[2015| l. 
Our ability to quantify the probability that a GRB occurred in a 
low-metallicity region of its globally metal-rich host galaxy will 
improve with larger sample sizes of nearby, high-metallicity, and 
preferably face-on host galaxies. However, for now, the observa¬ 
tions of the face-on host galaxy of GRB 020819B provide some 
support to the notion that GRBs can arise within metal-rich en¬ 
vironments. 

The increase in the detection rate of more dust and metal-rich 
GRB sightlines and host galaxies have already placed challenges 
on the standard collapsar model (e.g. Levesque et al.[ [2010bt 


to the host galaxy observations presented here. 

A BSG progenitor was also proposed in the case of 
GRB 111 209A ([Gendre et al.[ [2^ [Levan et al.[ [20T4l al¬ 
though see [Greiner et al. 2015 1 , which similar to GRB 101225A, 
had a very blue and compact host galaxy. There is cutTently 
limited information published on the host galaxy properties of 
GRB 121027A. Erom HST observations, [Levan et al.[ ( |2014[ ) es¬ 
timate the host galaxies of GRB 101225A and GRB 111209A to 
have radii of < 600 pc and ~ 700 pc respectively, and spectro¬ 
scopic observations of the host galaxy of GRB 111209A show 
that it had a sub-solar metallicity in the range 12H-log(0/H)=7.9- 
8.3 (Levan et al. 2014[ Rriihler et al.[ 2015[ l. The host galaxy 
of GRB 130925A, on the other hand, is much larger, with an 
effective radius of ~ 2.4 kpc, and a super solar metallicity of 
12H-log(0/H)- 8.8. 

In both the case of a WR star and BSG progenitor, the hy¬ 
drogen envelope has to be shed prior to the core-collapse with¬ 
out losing too much angular momentum. Within a more simple 
treatment of stellar evolution, where, for example, the rotation is 
modelled as a solid body, this places an upper limit on the pro¬ 
genitor metallicity of Z ^ 0.3 Zq (Hirschi et al. 2005 Woosley 
& Heger 2006[), which limits the mass-loss induced through 


Kriihler et al. [201 1| Savaglio et al. 2012| Perley et al. 201 3| l. 
However, in addition to the spectroscopically resolved observa¬ 
tions presented here, this host galaxy is distinctive in that it is the 
first example of a super-solar metallicity galaxy to host a GRB 
that makes up part of the cutTently small sample of ultra-long 
GRBs. As is the case for long-duration GRBs, the favoured the¬ 
oretical models for the ultra-long class of GRBs place an upper 
limit on the progenitor metallicity. 

5.2. Ultra-long GRBs; a distinct class of GRB? 

The prompt emission of GRB 130925A continued for an extraor¬ 
dinarily long time (~ 20 ks; [Evans et al.||20T4j[Piro et al.[[2014j l, 
and the X-ray afterglow light curve exhibited numerous large 


radiation pressure. Once the progenitor star increases to above 
^ 0.5^0, there is an excessive removal of angular momentum 
within the line-pressure driven stellar wind, and the blackhole- 
accretion disk system necessary to produce a GRB cannot form. 
However, in more complex models, additional considerations 
such as the effect of differential rotation or a companion star will 
also have an influence on the loss of angular momentum, and 
may thus relax the metallicity constraints. Eor example, when 
differential rotation is considered, it becomes possible for solar 
metallicity stars with initial stellar masses between 40 Mq and 
60 Mq to collapse to form a GRB (Georgy et al. |2012| l. 

A star may also retain a high angular momentum when it 
is part of a binary system either through tidal locking in which 
accretion onto the secondary star can significantly spin-up the 
stellar core, or through stellar core coalescence (jPodsiadlowski 


[et al.[[2004) l. Nevertheless, binary evolution models suggest that 
both cases are more likely to arise in low-metallicity environ¬ 
ments, and it is still not clear whether a progenitor core could 
retain enough angular momentum to form a GRB. Although the 
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core would be spun-up significantly during the helium burning 
phase, it is likely that a large fraction of this angular momentum 
would be lost prior to core-collapse. The gain in angular mo¬ 
mentum would result in an increased stellar wind, which in turn 
would reduce the radius of the star and increase the binary orbit, 
thus resulting in a weaker tidal interaction (Detmers et al. 2008| l. 
One possible solution is that the mass loss occurs preferentially 
along the rotational access, thus removing less angular momen¬ 
tum (Meynet & Maeder 2007| l. A weak coupling between the 
core and envelope of the progenitor star would also alleviate this 
constraint, such as in the absence of magnetic fields pzzard et al. 


|2004[|Petrovic et ^|2005[ [Detmers et al.||200^ . However, these 
additional conditions needed to form a GRB in a solar metallic- 
ity environment add similar complexities to the binary system as 
in the single star collapsar model. 

Owing to the very recent discovery of ultra-long GRBs, there 
has been little opportunity to investigate and develop possible 
progenitor models. Furthermore, given that GRB 130925A is the 
first example of an ultra-long GRB within a super-solar metal- 
licity host galaxy, little investigation has been done on how such 
ultra-long duration GRBs may form within solar metallicity en¬ 
vironments. The additional mass required to power an ultra-long 
GRB would imply that the conditions needed to maintain suffi¬ 
cient angular momentum at the time of core-collapse are likely 
to be even more stringent than those outlined above for standard 
long duration GRBs. 

An alternative, model-independent reconciliation to the ap¬ 
parent confiict between the low metal abundances predicted by 
collapsar models, and observations of host galaxies with solar 
metallicity and above, is that the high metallicities measured 
in the environments of some GRBs (either from afterglow ab¬ 
sorption lines or host galaxy emission lines) are not representa¬ 
tive of the progenitor metallicity, and the same argument could 
be used in the case of the ultra-long GRB 130925A. Although 
we measure a super-solar metallicity both along the disk of the 
host galaxy, which is consistent with the afterglow position of 
GRB 130925A, and in the outer regions, at SF2, the host galaxy 
of GRB 130925A is observed edge on, and we therefore cannot 
rule out a strong variation in metallicity along the radial direc¬ 
tion. If the GRB occurred at the far side of the host galaxy, it 
would be feasible for the GRB to have arisen in a metal-poor 
pocket of the host galaxy. Furthermore, HST observations of the 
host galaxy show evidence for an irregular galaxy morphology 
with a possible polar-ring structure ( |Tanvir et al.| |2013| l. Sig¬ 
nificant variation in the metallicity distribution may thus be ex¬ 
pected. Nevertheless, the detection of the very dusty and metal- 
rich host galaxy of GRB 130925A so early on after the discovery 
of the first ultra-long GRB would suggest that such environmen¬ 
tal conditions are not rare. 

Although there is substantial evidence indicating that long 
GRBs (and possibly also the ultra-long class) are more likely 
to arise in sub-solar metallicity environments with z < 0.5Zq 
(especially noticeable at z ^ 1.5; e.g. Modjaz et al. 2008 


Levesque et al. 2010a Graham & Frucn^f}' zul3 'Kmhlere 


al.| |2015|l, it is clear that long GRBs do not adhere to a strict 


metallicity cut-off. It is, therefore, important to consider the ef 
feet that other properties apart from metallicity have on the loss 
of angular momentum over a stars lifetime ([Levesque et al.[ 
2010c|l, such as anisotropic mass-loss (|Meynet & Maed^ 2007|l, 


wind clumping (Crowther et al. 2002| l, and magnetic processes 
( |Dessart et al.||2008| l. In order to compare observations with pro¬ 
genitor model predictions, it would be beneficial to quantify the 
probability that any particular progenitor model has in evolv¬ 
ing to form a GRB (e.g. |Hirschi et al.||2005[[van den Heuvel &| 
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[Portegies Zwart| |2013[ ), together with the observable signatures 
of any given model, wherever possible. This would then permit a 
direct comparison between the probability distributions and e.g. 
the observed metallicity distribution. However, ultimately, larger 
samples of spatially resolved GRB host galaxies are needed, and 
especially for host galaxies with super-solar metallicity, so that 
variations in the environmental conditions within the vicinity of 
the GRB relative to the rest of the host galaxy can be assessed, 
and underlying common properties can be identified. 


6. Summary 

In this paper we have presented a detailed analysis on the spec¬ 
troscopic X-shooter observations taken of the host galaxy of the 
ultra-long duration GRB 130925A (lasting for ^ 20 ks). The 
position of the slit fell along the disk of the galaxy, and fortu¬ 
itously crossed the GRB explosion sight, close to the nucleus of 
the host galaxy, and a star forming region located in the outskirts 
of the galaxy. Our observations suggest that the GRB occurred 
in a metal-rich region of the host galaxy, with a higher dust con¬ 
tent than that of a neighbouring HII region, which is in apparent 
conflict with the favoured progenitor models of both standard 
long duration GRBs, and the more recently discovered but more 
poorly understood class of ultra-long GRBs. 

Although the sample of long GRBs with dust-rich and super¬ 
solar metallicity host galaxies has drastically increased over the 
past half a decade, our spectroscopic data of GRB 130925A are 
only the second such example of spatially resolved observations 
of a super-solar GRB host galaxy, and they are the first such 
observations for an ultra-long GRB. The edge-on orientation of 
the host galaxy of GRB 130925A, and the very central sightline 
to the GRB through the galaxy disk, raises the possibility that 
GRB 130925A occurred within a metal-poor region, on the far 
side of the host galaxy. Nevertheless, the increased range in ob¬ 
served host galaxy properties provided by observations such as 
those presented here, and especially at such a relatively low red- 
shift, provides greater constraint on the GRB progenitor models, 
and ultimately, on the environmental properties traced by GRBs. 

With the currently small sample of ultra-long GRBs, the 
statistics for this class of event remain too small to draw any 
strong conclusions on the likely progenitors producing these 
incredibly long-lived and energetic events. Nevertheless, ob¬ 
servations such as our X-shooter data of the host galaxy of 
GRB 130925A, present new challenges to the favoured models. 
To continue in this vein it is thus important to be mindful of 
the dominant selection effects present in GRB and host galaxy 
samples, of which dust-extinction is possibly the most dominant, 
and to pursue observing programmes that address these selection 
effects, thus enabling us to quantify the extent of our biases. 
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